Effect of band filling in the Kondo lattice: A mean-field approach 
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The usual Kondo-lattice, including an antiferromagnetic exchange interaction between nearest- 
neighboring localized spins, is treated here in a mean-field scheme that introduces two mean-field 
parameters: one associated with the local Kondo effect, and the other related to the magnetic 
correlations between localized spins. Phases with short-range magnetic correlations or coexistence 
between those and the Kondo effect are obtained. By varying the number of electrons in the 
conduction band, we notice that the Kondo effect tends to be suppressed away from half filling, 
while magnetic correlations can survive if the Heisenberg coupling is strong enough. An enhanced 
linear coefficient of the specific heat is obtained at low temperatures in the metallic state. 
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I. INTRODUCTION 



The denomination Kondo latticea applies to systems 
which are characterized by a lattice of localized mag- 
netic moments coexisting with a conduction band. Ex- 
amples are found among Cerium compounds, such as 
CeAl 3 , CeCu6, CeCu2Si2, etc.,cl as .well as some Ura- 
nium or other rare-earth compounds .El At high tempera- 
tures, the localized moments behave essentially as inde- 
pendent impurities, similarly to what happens in dilute 
alloys. At low temperatures, a coherent heavy-fermion 
behavior is observed where the system resembles a Fermi- 
liquid with enhanced values of parameters such as the 
specific-heat constant 7 and the magnetic susceptibility 
X- This is a scenario found in the Cerium compounds 
mentioned above. However, deviations from this behav- 
ior are observed: some systems can show a magnetically 
ordered ground state, as happens, for example, in CeAl2, 
CeB 6 , and CePd2Al 3 , or become superconducting at low 
temperatures, as is the case of CeCu2Si2,ciQ CeCu2Ge2 
at high pressures ,□ and some .Uranium compounds like 
UPt 3 J2 UBei 3 J3 and URu 2 Si2,H among pQthers. Besides 
that, neutron-scattering experiments^! E3 point to the 
presence of strong short-range magnetic correlations in 
CeCu6, CeInCu2, and CeRu2Si2, for which a Fermi- liquid 
ground state is stable, but that are probably close to the 
conditions for a magnetic instability. 

This coexistence of magnetic correlations and Kondo 
effect w.as recently discussed by Iglesias, Lacroix, and 
Coqblir£3 within a mean-field approach that has points 
in common r-with a slave-boson treatment by Coleman 
and Andrei. Eil Their analysis was restricted to the case 
of a half-filled conduction band, which is of relevance for 
Kondo insulators $3 The role played by the conduction- 
band filling is obviously very important for metallic com- 
pounds that show heavy-fermion behavior, and to the 
discussion of underscreened systemsll3 

In this paper we study the effect of band filling in the 
stability of the Kondo state and short-range magnetic 
correlations in the Kondo lattice using a mean-field ap- 
proach closely related to the one employed in Ref. [L3. 



Two mean-field parameters are introduced, in connec- 
tion with the local correlations generated by the Kondo 
effect, and the non-local ones that indicate tendency to- 
wards magnetic order. These quantities behave as or- 
der parameters of thermodynamic phases where one of 
them or both are non-zero, while both vanish in the high- 
temperature "normal" phase. Besides the order parame- 
ters themselves, we calculate the free energy, the specific 
heat, and the single-particle density of states. We focus 
on the changes observed when one moves away from half- 
filling, as well as the competition between correlations 
induced by the Kondo coupling and those produced by a 
Heisenberg exchange term added to the original Kondo- 
lattice model. When this interaction is strong, the system 
first develops magnetic correlations, and only at a lower 
critical temperature does the Kondo effect appear, de- 
pleting but not suppressing the existing correlations. In 
contrast, when the Heisenberg interaction is small com- 
pared to the local Kondo coupling, the latter dominates, 
giving rise to a single phase where magnetic correlations 
are mostly induced by the Kondo effect, and can even 
change sign when the band filling is low enough. The 
specific heat that we obtain in the metallic case shows 
an enhanced linear coefficient, indicating that important 
qualitative aspects of the physics of heavy-fermion sys- 
tems are retained in the mean-field approach. 

In Sec. we introduce the model Hamiltonian and 
the fermionic representation of localized moments. In 
Sec. [Hi] we choose the relevant fields in terms of which 
we rewrite the Hamiltonian, perform the mean-field de- 
coupling, find the energy eigenvalues, and introduce self- 
consistency equations for the mean-field parameters, as 
well as relations to obtain the physical qu antities of in- 
terest. Our results are presented in Sec. |IV[ and a critical 
discussion of them appears in Sec. [Vj 



II. MODEL HAMILTONIAN 



The usual Kondo-lattice Hamiltonian is 
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H 
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(1) 



where s i stands for the total conduction-electron spin at 
lattice site i, S i is the spin operator associated to the 
localized moments, and the first term in the right-hand 
side describes the conduction band in usual notation. We 
choose spin-1/2 localized moments, assigning them to "/ 
electrons". We also add a Heisenberg-like interaction 
between nearest-neighbor localized spins. The Kondo- 
lattice Hamiltonian then reads 
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(2) 



The term with the factor Eq is just a constant provided 
that we remain in the subspace of unit occupation num- 
ber for the / electrons at every site, which is the con- 
straint that ensures the equivalence of Eqs. ([!]) and (0) if 
the exchange interaction is also included in the original 
model. 

The spin operators are rewritten in the usual fermionic 
representation 
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There is no unique way of rewriting the Hamiltonian in 
terms of fermion operators. For instance, the terms in- 
volving the z component of the spins can be left as in 
Eq. (0), or the number operators can be explicitly writ- 
ten as products of creation and annihilation operators, 
whose ordering can then be altered using the fermionic 
anticommutation relations. The final choice of form to 
write down these terms will actually depend on the mean- 
field decoupling that will be employed, as we will discuss 
in the next section. 



III. MEAN-FIELD SCHEME 

The first step in constructing the mean-field Hamilto- 
nian is to choose the relevant fields. The original Hamil- 
tonian is naturally written in terms of spin fields. How- 
ever, if we just introduced a mean-field decoupling of spin 
products we would end up with a magnetic Hamiltonian. 
Since we intend to describe the Kondo effect as well, we 
have to choose fields that couple c and / electrons at 
the same site, since these will tend to form local singlets. 
We can also couple / electrons from neighboring sites in 
order to describe the short-range magnetic correlations 
that will develop between the localized moments. A con- 
venient way to do this is to define the (Hermitian) fields 
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where a labels the spin orientation (up or down) with re- 
spect to the z axis, and i and j are nearest-neighbor sites 
in the definition of r, ijCT . It is straightforward to write 
down the contributions of transversal spin components 
in terms of these fields. We have 



S\ Sj + Sf Sj = — f ij a t" 



(5) 



where we have made explicit usage of the constraint 
n{ = 1. Products involving the z component of the spins 
are slightly more complicated to deal with, and some- 
times these are decoupled as spin fields.Eij However, this 
breaks the spin symmetry, and results in the stability of 
a magnetic state, with total suppression of the Kondo ef- 
fect, unless the local average magnetic moment is forced 
to be zero. Hence, we will avoid this procedure, using 
the same kind of representation for all spin components. 
After some algebraic manipulations, we obtain 
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Due to the constraint n\ = 1 at all sites, the terms con- 
taining number operators in Eqs. (|^) will yield a con- 
stant shift of the energies (after summation over the lat- 
tice sites). This terms can then be dropped, leaving the 
Kondo and Heiscnbcrg parts of the Hamiltonian respec- 
tively written as 
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and 
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(8) 



We now proceed with a standard mean-field decoupling 
of the operator products in these Hamiltonians, introduc- 
ing the mean-fields Xi a = (Xi a ) and r t3CT = (Fy CT )' From 
Eqs. (Q) and (||) ones sees that each field couples with a 
spin-independent quantity. Thus, there will be no break- 
down of spin symmetry, i. e., no magnetic states. Taking 
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this into account, together with translational invariance, 
we can simplify the notation, using A^ = A and r,y CT = T. 

We are now in a position to put all contributions to- 
gether: the noninteracting terms of Eq. (^|) and the mean- 
field versions of the interaction Hamiltonians (Q) and (j^). 
Reverting to the fermionic representation, and adding a 
chemical potential /i for the conduction electrons, we can 
write the complete mean-field Hamiltonian 
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where z represents the coordination number of a site in 
the lattice, and N is the total number of sites. 

We will choose the conduction band to be a tight- 
binding band of width 2W (with nearest-neighbor hop- 
ping only) in a (hyper)cubic lattice in d dimensions, so 
that we can write 
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where [i labels the wavevector components, and a is the 
lattice parameter. With this choice, the term with Jh in 
Eq. (^) will have the same k-dependence as the conduc- 
tion band, except for the band width. We see, thus, that 
the mean-field decoupling introduces a non-zero band 
width to the / electrons. At this point, we have to re- 
place the constraint n\ = 1 with the much weaker con- 
straint (n{) = 1. Hence, the virtual charge fluctuations 
introduced by the fermion representation of the localized 
spins have been "put in the mass shell" by the mean- 
field decoupling. Of course, this real charge fluctuations 
are spurious, but the method relies on the expectation 
that the physics of the system will be reasonably pre- 
served on the phase-space surface where the constraint 
( n i) ~ 1 holds. In order to improve on this point one 
would have to envisage better ways to take into account 
the constraint of single occupancy on the / levels at each 
site. 

Leaving aside, for the moment, the constant terms i?x 



and En, the Hamiltonian (E 
ing the energy eigenvalucsE 



is easily diagonalized, yield- 
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where B = —AdJjjT/W. These energies are spin- 
degenerate, as we remarked before. With them and the 
corresponding eigenvectors we can compute any average 
of relevant physical quantities in the system. The chem- 
ical potential /i is determined by the equality (n k ) = n, 
where n is the (chosen) density of conduction electrons. 
Eq plays the role of a chemical potential for the / elec- 
trons, being fixed by the condition (ri£) = 1. The mean- 
field parameters A and T are determined through the self- 
consistency equations 
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We also check the self-consistency process by evaluat- 
ing the free-energy, which can be written as 



F = -2T ^ log l + e- B £ /T 

k,a=± 

+ E K + E H - (E - fin)N 



(14) 



where T is the temperature (in energy units). Finally, 
we can calculate the average internal energy 

E = 2 EZf(EZ)+E K + E u -(E -tm)N , (15) 

k,a=± 

where /(e) stands for the Fermi function. The last term 
in Eqs. (|l4) and ( |l5| ) is needed to compensate for the fact 
that fi and Eq are already included in the Hamiltonian. 
The specific heat can now be obtained as c v = OE/dT, 
and we can check for an enhanced linear coefficient 7, 
which is a signature of heavy-fermion behavior. 



IV. RESULTS 

In order to be close to real systems, we considered a 
three-dimensional simple-cubic lattice (z = 6). However, 
given that all wave- vector dependence occurs through £k, 
we turn all k sums into integrals over the conduction- 
band energies, for which we use a semi-elliptical density 
of states of the form 



V{e) 



1(1 
0. 



-1 < £< 1 

otherwise . 



(16) 



Here we have set W — 1 , so that all energies are measured 
in units of the half band width. We choose both Jk and 
Jh to be negative, which corresponds to the usual Kondo 
coupling, and antiferromagnetic Heisenberg interactions 
between nearest-neighboring localized spins. The latter 
would be naturally interpreted as generated by a superex- 
change mechanism, and could also mimic the RKKY in- 
teraction in the mean-field Hamiltonian, although this is 
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a debatable issue: we will see below that there exist in- 
duced spin correlations due to the Kondo coupling alone, 
and these are more closely related to the RKKY mecha- 
nism. 
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FIG. 2. Same as in Fig. |l| but with a large Heisenberg 
interaction: J K /W = -0.5, J K /W = -1.0. 



FIG. 1. Mean-field parameters A and F as functions of tem- 
perature for n — 1 (top) and n = 0.6 (bottom), in the case 
where the Heisenberg interaction is small compared to the 
Kondo coupling: J K /W = -0.5, Jn/W = -0.1. 

After numerically solving the set of self-consistency 
equations, we find that both A and T behave as typical 
mean-field order parameters, which vanish at high tem- 
peratures, and become non-zero through second-order 
phase transitions at well defined critical temperatures. 
We identify the region where A / with the Kondo 
regime, and the corresponding critical temperature is in- 
terpreted as the Kondo temperature (Tk) for the lat- 
tice. We call Tborr the temperature below which short- 
range magnetic correlations appear (r ^ 0). When 
|Jh| << \Jk\, both transitions occur at the same tem- 
perature, while |«7h| comparable to or larger than |Jk| 
giveSj-Jcorr > Tk, as obtained earlier for the half-filled 
case.tj We show a typical example of this behavior in 
Figs. |l| and [| where we compare the results for a half- 
filled and a less-than-half-filled conduction band. (Due 
to particle-hole symmetry, the results are actually sym- 
metric around n = 1.) 



Departure from the half-filled situation causes a re- 
duction of the critical temperatures, except for Tc rr in 
the large- 1 Jh | regime. Although the sign of T does not 
have a direct physical meaning, since spin correlations 
are related to the square of the operator associated to T 
according to Eq. (6|), we believe that the change of sign 
observed in Fig. [l for the case of low band filling indi- 
cates that the nearest-neighboring spin correlations have 
become ferromagnetic, as expected in this limit. This 
is only observed at low temperatures, where the Kondo 
effect is large, and the induced part of the correlations 
dominate. In contrast, close to Tk, A is small, and the an- 
tiferromagnetic correlations due to the Heisenberg inter- 
action dominate. In order to check this, we would have to 
calculate the average (rfj CT ) that appears in the spin-spin 
correlations, Eq. (|g). Unfortunately, this is related to av- 
erages of products of number operators in each site, which 
are non-trivial in spite of the non-interacting nature of 
the effective Hamiltonian. More specifically, these prod- 
ucts involve number operators at displaced wavevectors, 
which prevents the transformation of momentum summa- 
tions into energy integrals. Notice, however, that when 
r goes to zero, spins in different sites become decoupled 
in the mean-field Hamiltonian. Then, their correlations 
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go to zero, which should happen when they change from 
antiferro- to ferromagnetic. 
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FIG. 3. Variation of the mean-field parameters with band 
filling at fixed temperature (T/W = 0.05). In the left panel 
we have dominant magnetic correlations (Jk/W = —0.5, 
Jh/W = —1.0), while in the right panel the Kondo effect 
dominates {Jk/W = -0.5, Jh/W = -0.1). 
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tions of the band filling for a fixed (low) temperature 
(T/W=0.05). For large \ J H \ (left panel), the Heisenberg 
coupling dominates, keeping T almost unchanged, while A 
is strongly reduced as we move away from n = 1, since we 
have an underscreening situation. This picture is com- 
pletely changed in the small- |Jh| regime (right panel). 
Now both parameters are reduced as one moves away 
from half-filling. One sees that T is now equally sensi- 
tive to n, which suggests that in this regime the part 
of the inter-site correlations induced by the conduction 
electrons is dominant. It eventually yields ferromagnetic 
spin correlations, as discussed above. We checked that a 
similar behavior is observed if we fix Jh = 0, when all 
spin correlations are induced by the conduction electrons. 

We have cut the curves in the right panel of Fig. || for 
small n because the convergence of our numerical calcu- 
lations, which involve iterating the set of self-consistency 
equations for the mean-field parameters and occupation 
numbers, becomes very delicate in that region. In order 
to see what is happening, we evaluate the free energy F 
[Eq. (|l4])] for fixed A, T, n, and T, adjusting only \i and 
Eq. Typical contour plots of F at low temperature as a 
function of A and T, for half-filling and for n very small, 
in the small- |Jh| regime, are shown in Fig. |]. We can see 
that the minimum at finite A and T that exists at half- 
filling is displaced towards small values of A and negative 
values of T as n is reduced. Eventually, a new minimum 
appears for T > 0, since the energy eigenvalues, Eq. (p^), 
are symmetric with respect to the sign of T when A = 0. 
Thus, when A approaches zero (as n is reduced) we have 
two shallow minima separated by a very broad and low 
maximum at A = T = 0, which causes problems to the 
numerical iteration of the self-consistency equations since 
distant points in the A-r space are very close in free- 
energy. 





FIG. 4. Equal free-energy contours at n = 1 (top) and 
n = 0.1 (bottom) for | Jh| small, showing the complex topolgy 
that arises at low filling. The values of the free energy for 
each contour, quoted on the right, are in units of the half 
band width. 

In Fig. p| we plot the mean-field parameters as func- 



FIG. 5. Variation of the specific heat (per unit cell) with 
temperature for n = 1 (continuous) and n = 0.6 (dashed), 
when the Heisenberg interaction is weak (Jk/W = —0.5, 
Jh/W = —0.1). Notice the insulating and metallic behav- 
ior at low temperatures in each case, and the lambda-shaped 
transition at the point where both the Kondo effect and mag- 
netic correlations disappear. 
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FIG. 6. Same as in Fig. pi but for strong Heisenberg cou- 
pling (Jk/W — —0.5, Jh/W — —1.0), showing the transitions 
at different temperatures. 

Our results for the specific heat as a function of tem- 
perature are presented in Figs. |^ and ^ The phase tran- 
sitions appear in the usual "lambda shape" characteristic 
of mean-field approximations. There are two jumps in the 
case of large | Jh | , at the temperatures where the Kondo 
effect and the magnetic correlations disappear. These 
points coincide when | Jh| is small. We can also see the 
existence of a gap at the Fermi level for the half-filled 
band, yielding an exponentially vanishing specific heat 
for T — > 0. In contrast, away from half- filling one sees 
a linear regime at low temperatures. The linear coeffi- 
cient obtained for n = 0.6, \Jk\ — 0.5, and |Jh| = 0.01 is 
7 ~ IO70, where 70 is the corresponding value in the ab- 
sence of interactions. We have, thus, some enhancement 
of 7, although the true heavy-fermion limit (7 ~ 10 2 7o) 
is not reached. Such an enhancement should correspond 
to a high density of states (DOS) at the Fermi level. We 
calculated the DOS as 



p(w) 



in terms of the one-particle Green's functions 



1 



z-Ei 



(17) 



(18) 



The DOS depends on temperature through A and T, 
which appear in the energies Eg. We plot the low- 
temperature DOS in Fig. |?], for both n — 1 and n = 0.6, 
in the small- |Jh| regime. We can see that the Fermi level 
(u; = 0) falls in the gap for n = 1, while it is located at a 
point corresponding to a very high density of states away 
from half filling. The gap that is observed immediately 
above the Fermi level in the latter case explains the re- 
duction (or at least levelling off) of the specific heat just 
above the linear region in Fig. H. As the temperature is 
increased, the gap in the DOS as well as the peaks around 
it are narrowed until a superposition of a free conduction 
band and a localized /-level is recovered above Tk. 
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FIG. 7. Total single-particle density of states at 
T/W = 0.01, for n = 1 (top) and n = 0.6 (bottom), with 
Jk/W = —0.5, Jh/W — —0.1. The frequencies are measured 
in units of W. Notice that the Fermi level (vertical dashed 
line) lies in the gap at half-filling, and inside a high peak for 
n < 1. 



V. CONCLUSIONS 



We have presented here a detailed mean-field anal- 
ysis of the Kondo lattice, emphasizing the effects of 
conduction-band filling in the properties of the system. 
The regime where the Kondo effect is observed, and 
the one in which short-range magnetic correlations are 
present appear as thermodynamic phases with well de- 
fined critical temperatures. These phases are character- 
ized by non-zero values of the mean-field parameters as- 
sociated to local correlations between localized and con- 
duction electrons (A), and between localized electrons 
on nearest-neighboring sites (r). We find a regime in 
which the Kondo effect and short-range magnetic cor- 
relations coexist, in agrrement with—experimental ob- 
servations in some Ce compounds.Bilj The main con- 
sequence of changing the electronic density in the con- 
duction band with respect to half-filling is a tendency 
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to suppress the Kondo effect. Magnetic correlations are 
equally suppressed if a direct Heiscnberg interaction be- 
tween nearest-neighboring localized spins is weak, while 
they are almost insensitive to the band filling for strong 
Heisenberg interactions. In the latter case the critical 
temperatures for the two order parameters are different, 
and one notices a partial suppression of magnetic corre- 
lations in the region where the Kondo effect is observed. 
The phase transitions appear as lambda-shaped peaks in 
the specific heat, which also shows an enhanced linear co- 
efficient in the metallic situation (n/ 1), corresponding 
to an enhanced density of states at the Fermi level. 

To a good extent, the physics of Kondo systems is re- 
covered here, at least qualitatively, both for insulators 
(n = 1) and conductors (n / 1). A weak point of this 
treatment is the impossibility to describe magnetically 
ordered states. Here one would have to go back to the 
choice of relevant fields, trying to write down a Hamil- 
tonian in which both Kondo and spin fields were present 
before the mean-field decoupling. It would be interesting 
to improve the method along these lines without loosing 
its formal simplicity. This simplicity constitutes the main 
quality of this approach, allowing evaluation of many rel- 
evant physical quantities in an unrestricted range of pa- 
rameters and conditions, and yielding qualitative under- 
standing of the behavior of the system as a whole. 
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